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ABSTRACT 


Epichloé species are fungal symbionts (endophytes) of grasses, many of which are benign or mutualistic and have a 
balance of horizontal (contagious) and vertical (seed-borne) transmission, whereas others mainly transmit horizontally 
and are more antagonistic. Over the past eight years several Epichloé species have been described based largely on 
the biological species concept. We conducted a multi-gene phylogenetic analysis to evaluate these endophytes as 
phylogenetic species, and thereby assess the relationship of phylogenetic and biological species. Variation mainly in 
introns of genes encoding B-tubulin (tub2), translation elongation factor l-a (tefl), and actin (act1) provided robust 
phylogenetic signal distinguishing the described Epichloé species. Outgroup rooting split the genus into two major 
groups. One group included most species with balanced transmission strategy, and in this group the phylogenetic and 
biological species concepts corresponded well. In contrast, these species concepts poorly corresponded for the other 
group, the Epichloé typhina complex, with predominantly antagonistic, horizontally transmitted endophytes. We suggest 
that the balance of vertical and horizontal transmission may promote ecological (host) specialization and subsequent 
genetic isolation as mechanisms promoting speciation; whereas strict horizontal transmission may select for broader 
host ranges, slow the development of genetically isolated species, and thereby increase lineage sorting effects that 
cause conflicts between phylogenetic and biological species. 

Key words: biological species, Clavicipitaceae, Epichloe, grass endophytes, molecular phylogenetics, phylogenetic 
species, Poaceae, symbiosis. 
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Fungi in the family Clavicipitaceae (Ascomycota) 
are a diverse group of organisms capable of inter- 
acting with an even broader assemblage of hosts. 
Contained within the Clavicipitaceae are fungal 
and invertebrate pathogens (tribe Cordycipeae), 
and plant pathogens (tribes Ustilaginoideae, Clav- 
icipeae, and Balansieae) (Diehl, 1950). Members of 
the tropical tribe Ustilaginoideae are biotrophs of 
bamboos, forming hardened sclerotic stromata on 
stems (Diehl, 1950), while the Clavicipeae destruc- 
tively infect grass florets and replace host ovules 
with fungal sclerotia (Tudzynski et al., 1995). Fungi 
in the tribe Balansieae form perennial infections of 
grasses, yet cause little or no symptoms of disease. 
The extent of host colonization by members of this 
tribe varies from fairly localized, epibiotic infec- 
tions to a more systemic, endophytic growth habit 
(Glenn et al., 1996). Species in the genus Epichloé 
lie at this endophytic extreme, systemically infect- 
ing cool-season grasses in the subfamily Pooideae 
(Schardl & Phillips, 1997). All Epichloé species are 


intercellular inhabitants of their host grasses, caus- 


ing no visible symptoms of infection for the majority 
of the symbiont and host life cycles (Fig. 1). The 
only macroscopic sign of infection occurs as the 
grass reaches reproductive maturity and initiates 
flowering. At this stage, the endophytes may cause 
abortion of some or all inflorescences (grass-choke 
disease), and replace them with their own sexual 
structures (stromata). When a stroma is fertilized 
by spermatia from a stroma of an opposite mating 
type, the sexual state develops and culminates in 
production and ejection of contagious spores. The 
level of antagonism to the host is directly related 
to the level of fungal sexual expression, as a "strong 
choker" effectively shuts down host seed produc- 
tion (Schardl et al., 1997). While some Epichloë 
species seem to choke (stromatize) every develop- 
ing host inflorescence, other species reproduce sex- 
ually on only a portion of the host tillers, the other 
inflorescences develop normally, and the endophyte 
enters into the developing ovule of the host seed. 
The result is vertical (asexual) transmission of the 
endophyte from the maternal host to her offspring. 
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Figure 1. Generalized life cycle of Epichloé species using E. festucae іп Festuca rubra as an example. Both trans- 
mission life cycles are exhibited by several Epichloé species (see text). Some Epichloé species (such as E. typhina) 
exhibit little or no asexual transmission in host seeds and are, therefore, predominantly transmitted horizontally via 
ascospores. Asexual derivatives of Epichloé species, conventionally classified in genus Neotyphodium, rely entirely on 
vertical transmission for their dissemination. Modified from Schardl et al. (1997). 


Previous research based mainly on artificial mat- 
ing tests has capitalized on the presence of the Ep- 
ichloé sexual cycle and heterothallism (obligate out- 
crossing; White & Bultman, 1987) to circumscribe 
interfertility groups (mating populations) that cor- 
respond to, or approximate, biological species 
(Leuchtmann et al., 1994). Here we use, as a work- 
ing definition of biological species, populations 
whose members have the potential to interbreed in 
nature to produce viable offspring, but who cannot 
interbreed with members of other species (Claridge 
et al., 1997; Mayr, 1940). To date, nine mating pop- 
ulations of Epichloé have been identified (Schardl 
et al., 1997). 

Many grass endophytes appear to be restricted 
to an asexual lifestyle, and are only transmitted ver- 
tically (Sampson, 1933). In fact, the presence of 
such endophytes was widely recognized only after 
the realization that forage grasses infected with 
these asexual species were toxic to livestock such 
as cattle and sheep (Bacon et al., 1977). Subse- 
quent studies have shown these endophytes to be 
responsible for producing novel secondary metab- 
olites responsible for the toxicoses (Kennedy & 
Bush, 1983), as well as conferring on their hosts 
increased drought tolerance (Arechavaleta et al., 
1989) and many other dramatic fitness enhance- 
ments (Schardl & Phillips, 1997). Despite the ap- 
parent inability of these asexual endophytes to 
spread contagiously, they can be abundant because 


they are seed transmitted at nearly 100% efficiency 
(Siegel et al., 1984). Although their hyphal mor- 
phology, conidiophore size and structure, and host 
associations all indicate their relatedness to Epich- 
loé species, their lack of a sexual phase has pre- 
cluded their classification based on a biological 
species concept. 

The development of modern molecular tech- 
niques such as polymerase chain reaction (PCR) 
and DNA sequencing allows classification based on 
genetic relationships (phylogenies) to be consid- 
ered. The strength of using gene sequence data for 
inferring relationships arises from a wealth of pre- 
cisely comparable characters (aligned positions of 
nucleotide sequences) between even the most dis- 
tantly related species (Felsenstein, 1988). Relation- 
ships emerging from such analyses have, in many 
cases, corroborated classification schemes based on 
morphological characteristics (Hillis, 1987). Molec- 
ular data have also been used to either strengthen 
or weaken support for hypotheses regarding organ- 
ismal and cellular evolution (Hillis, 1987). For ex- 
ample, gene sequence data indicate that cellular 
organelles such as mitochondria and chloroplasts 
are most closely related to some free-living bacte- 
rial species, lending credence to the endosymbiont 
theory of the modern-day eukaryote cell as a com- 
posite of two or more once distinct prokaryote spe- 
cies (Margulis, 1991). 

Gene sequence data also allows phylogenetic 
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inferences for organisms that lack distinguishing 
features for reliable morphological classification. 
Because of rather little morphological variation 
many sexual endophytes were originally classified 
as Epichloé typhina. Recently, a number of species 
have been described based on molecular phyloge- 
netic evidence, interfertility groups, and morphol- 
ogy (Leuchtmann & Schardl, 1998; Leuchtmann et 
al, 1994; Schardl & Leuchtmann, 1999; White, 
1993, 1994). A phylogeny has been inferred from 
partial sequences of 26S large subunit rRNA genes 
from Epichloé species and other members of the 
family Clavicipitaceae (Kuldau et al., 1997). Those 
endophytes assigned to the genus Epichloé through 
morphological data and mating tests all formed a 
monophyletic clade within the Clavicipitaceae. A 
more comprehensive analysis within the genus Ep- 
ichloé was conducted using non-coding regions of 
the B-tubulin gene (tub2) (Schardl et al., 1997). 
Each sexual isolate contained a single tub2 copy, 
reflecting the haploid state of the Epichloé myce- 
lium. Most but not all mating populations corre- 
spond to well defined clades in phylogenetic anal- 
yses based mainly on tub2 introns. However, more 
detailed analyses are needed to further test the in- 
ferred relationships. 

Similar analysis of tub2 sequence data obtained 
from the strictly asexual grass endophytes revealed 
quite novel findings. The majority of these isolates 
actually appeared to have multiple tub2 copies 
(Schardl et al., 1994; Tsai et al., 1994; Moon et al., 
2000; K.D.C. & C.L.S., unpublished obs.). Through 
a combination of restriction endonuclease digestion 
and selective priming in PCR, as well as by clon- 
ing, these copies could be isolated and sequenced 
independently. Inclusion of these sequences in a 
phylogenetic analysis of Epichloé species revealed 
that each tub2 copy contained in the genome of an 
asexual grass endophyte aligned most closely with 
that of a different sexual Epichloé species (Schardl 
et al., 1994; Tsai et al., 1994). For example, Neo- 
typhodium coenophialum, a common asexual en- 
dophyte of tall fescue [Lolium arundinaceum 
(Schreb.) S. J. Darbyshire — Festuca arundinacea 
Schreb.] was found to harbor three tub2 copies cor- 
responding with those of three different Epichloé 
species: E. festucae, E. typhina, and E. baconii 
(Tsai et al., .994). Similar results were obtained by 
analysis of ther genes (specifically, those utilized 
in this paper; K.D.C. & C.L.S., unpublished data). 
Furthermore, other asexual endophytes from Fes- 
tuca L. spp., Lolium L. spp., Poa L. spp., and Stipa 
L. spp. have genes derived from different combi- 
nations of Epichloé species, such as E. baconii with 


E. typhina (Tsai et al., 1994), E. baconii with E. 
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bromicola (Moon et al., 2000), E. amarillans with 
E. typhina (Schardl & Wilkinson, 2000), E. festucae 
with E. typhina (K.D.C., unpublished), E. festucae 
with E. elymi, E. bromicola with E. typhina, and E. 
bromicola with E. festucae (K.D.C., unpublished). 
Based on these findings, most asexual endophytes 
related to Epichloé spp. are likely to be genetic 
hybrids arising through the parasexual process of 
hyphal anastomosis (Schardl et al., 1994; Schardl 
& Wilkinson, 2000), and the asexual tall fescue 
endophyte is apparently derived from at least two 
such hybridizations. The genus Neotyphodium has 
been adopted to contain such endophytes with close 
relationships to Epichloé spp., but lacking a sexual 
stage in their life cycles (Glenn et al., 1996). A 
morphological difference that often distinguishes 
Neotyphodium spp. from their Epichloé relatives is 
that the asexual endophytes tend to have larger co- 
nidia, probably as a result of increased genome 
content due to hybrid origins (Kuldau et al., 1999), 

A potential problem that arises when sequence 
data are used to construct phylogenies is whether 
or not the gene tree accurately represents the true 
evolutionary relationships between species (species 
tree). For example, similarities between two se- 
quences can arise from convergence, gene conver- 
sion, or homology (common ancestry). Homology 
can further be divided into orthology (common an- 
cestry tracing back to a speciation event) and par- 
alogy (common ancestry tracing back to a gene du- 
plication event) (Moritz & Hillis, 1996). Sampling 
of paralogous sequences can cause incorrect infer- 
ences about species phylogeny, so only orthologous 
sequences should be evaluated. Therefore, the 
genes chosen in such analyses are of paramount 
importance if one wishes to obtain an accurate evo- 
lutionary picture. 

Another way of measuring the accuracy of a gene 
tree in reflecting the actual species tree is to sam- 
ple other genes. If trees developed from different, 
independently assorting genes contain identical 
branching orders and taxa placement, this can be 
regarded as strong support for these estimations as 
accurately revealing the evolutionary relationships 
among the species being studied. An additional ad- 
vantage to multi-gene phylogenetics is that popu- 
lations that are, or have been, interbreeding can be 
distinguished from those that are genetically iso- 
lated in nature (Geiser et al., 1998; Burt et al., 
1996). This is important for fungal evolutionary bi- 
ology because on the one hand many lack a known 
sexual stage, and on the other hand the ability to 
mate in experimental crosses (biological species) 
does not guarantee that the process is relevant to 
natural populations. This paper reports the char- 
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acterization and analysis of two more genes in the 
Epichloé genome. Portions of the translation elon- 
gation factor l-a (tefl) gene (O'Donnell et al., 
1998) and the actin (act1) gene (Fidel et al., 1988), 
with variable introns, were amplified and se- 
quenced for this group of closely related grass en- 
dophytes. The sequence data were used to con- 
struct separate gene phylogenies, which could be 
compared to those obtained from tub2 sequences. 
The results indicate different phylogenetic patterns 
for symbiont species that themselves differ in rel- 
ative mutualistic or antagonistic character. 


MATERIALS AND METHODS 
BIOLOGICAL MATERIALS 


Fungal isolates, their hosts, and sequence acces- 
sion numbers are listed in Table 1. Endophytes 
were isolated by plating infected, surface-sterilized 
leaf sheath and stem material onto potato dextrose 
agar (PDA, Difco Laboratories, Detroit, Michigan) 
containing streptomycin and penicillin (1 mg L*' 
each). Plates were checked regularly for contami- 
nation and endophytic growth. Most strains grew 
out 5—14 days after plating. Isolates were stored in 
50% glycerol at —80°C (Schardl & Ап, 1993) and 
on PDA slant cultures overlain with autoclaved 
light mineral oil. 


MATING TESTS 


Plants naturally infected with Epichloé were ver- 
nalized over winter at an experimental garden in 
Zurich, Switzerland, to induce stromata formation. 
The following spring, plants were moved into a 
greenhouse where mating tests were performed on 
emergent stromata. Spermatia were transferred be- 
tween mating types by rubbing stromata of opposite 
mating type together, or by applying cultures con- 
taining conidia onto stromata. Successful matings 
between opposite mating types (of the same mating 
population) were evidenced by the production of 
tan to orange perithecia after 3—6 weeks, from 
which ascospores could be collected. Single asco- 
spore isolates were grown on media and subjected 
to isozyme analysis. Independent assortment of pa- 
rental isozyme alleles in ascospore progeny con- 
firmed successful matings. 


PREPARATION OF GENOMIC DNA 


Subcultures were generated by grinding a small 
amount of fungal mycelium in sterile water using 
an autoclaved mini-pestle and spreading approxi- 
mately 300 pl of this suspension atop sterile cel- 
lophane disks on PDA plates. The permeable cel- 
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lophane allowed hyphal growth and facilitated 
subsequent harvest of mycelium. Cultures were in- 
cubated at room temperature for two weeks, after 
which mycelium was harvested and freeze dried. 
Total fungal DNA was prepared as previously de- 
scribed (Byrd et al., 1990; Schardl et al., 1997) or, 
preferably, the method of Al Samarrai and Schmid 
(2000). 


ASSESSMENT OF GENE LINKAGE 


To assess the possibility of linkage between the 
three genes used in this analysis, spermatia from iso- 
late ATCC200736 (Epichloé typhina from Lolium per- 
enne) were used to fertilize stromata of isolate 
ATCC200740 (E. typhina from Dactylis glomerata). 
Since these two E. typhina isolates are of opposite 
mating type, perithecia developed from which asco- 
spores were obtained. Inheritance of parental alleles 
in each progeny was assessed by polymerase chain 
reaction (PCR) with primers specific for each tub2, 
tefl, and асї1 allele. All primers are listed in Table 
2. Amplification of tub2 spanning introns 1-3 was 
performed using non-selective primer tub2-exon]d—1 
with selective primers tub2-intron3u-200736 (selec- 
tive for E. typhina ATCC200736) and tub2-intron3u- 
200740 (selective for E. typhina ATCC200740). The 
non-selective primer annealed to the exon 1-intron 1 
junction while both selective primers annealed to 
complementary regions within intron 3. A region of 
the act] gene approximately 300 bp long (including 
introns 2-4) was amplified using two selective primers 
for each parental isolate. Primers actl-intron2d- 
200736 and act] -intron-4u-200736, respectively, an- 
nealed to sites within introns 2 and 4 of the 
ATCC200736 act1 allele. Similarly, primers act1-in- 
tron2d-200740 and actl-intron4u-200740 annealed 
to these same regions of the ATCC200740 асі] allele. 
Partial amplification of the tefl gene, including in- 
trons 1—2, was performed using the non-specific prim- 
ег tefl-exonld with specific primers tefl-intron3u- 
200736 for ATCC200736 or tef1-intron3u-200740 for 
ATCC200740. The non-specific primer anneals to 
conserved sequence within exon 1, while both isolate 
specific primers anneal to complementary sequence 
within intron 4. PCR was performed in 50 pl reac- 
tions containing 0.2 uM of each deoxynucleotide tri- 
phosphate (dATP, dCTP, dGTP, and dTTP; Boehringer 
Mannheim, Indianapolis, Indiana), 1 X PCR buffer 
(10 mM tris(hydroxymethyl)aminomethane hydrochlo- 
ride [Tris-HCl] pH 8.3, 15 mM MgCl, 500 nM of 
each primer (see below), 10 ng total fungal DNA, and 
0.02 U pl"! of AmpliTaq gold DNA polymerase (Рег- 
kin-Elmer, Foster City, California). The thermal cycler 
(Perkin-Elmer model 2400) was programmed as fol- 
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Table 2. 


Oligonucleotide 
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Oligonucleotide primers used in this study. Primers are listed in order of appearance in text. 


designation Gene Sequence (5'—3') Position Orientation 
tub2-exonld-1 tub? GAGAAAATGCGTGAGATTGT exon l-intron l junction downstream 
tub2-intron3u-200736 tub2 = GCATGATAACACTTTTCCTTTCCA intron З upstream 
tub2-intron3u-200740 — tub2 ССАТСАТААСАССТТТТАТТТССТ intron З upstream 
actl-intron2d-200736 асг] GCGTAGACCCACGCA ACGTA intron 2 downstream 
actl-intron4u-200736 — act! GTGCATGTTTGCCAGCCAGA intron 4 upstream 
actl-intron2d-200740 асі] GCATAGACCCACGCA ACGTG intron 2 downstream 
actl-intron4u-200740 act] GTGCATGTTTGCCAGCCAGT intron 4 upstream 
tefl-exonld tefl GGGTAAGGACGAAAAGACTCA exon | downstream 
tefl-intron3u-200736 tefl ATCACGCGAATGTCATGTTG intron 4 upstream 
tefl-intron3u-200740 tefl ATCACGCGAATGTCATGTTA intron 4 upstream 
tub2-exon4u-1 tub2 ТССТСААССАССТСАССАСС ехоп 4 upstream 
tefl-exon6u-1 tefl CGGCAGCG ATA ATCAGGATAG exon 6 upstream 
actl-exonld-1 act! — TAATCAGTCACATGGAGGGT exonl-intron | junction downstream 
actl-exon6u-] actl AACCACCGATCCAGACAGAGT exon 6 upstream 
actl-exon3u-1 act] CTTGCACATGCCCGA ACTACC exon 3 upstream 
actl-exon3d-1 actl GGTAGTTCGGGCATGTGCAAG exon 3 downstream 





lows: one incubation at 94°C for 9 min.; then 40 cy- 
cles of 94°C for 60 sec., 60°C for 60 sec., and 72°C 
for 60 sec. 


SEQUENCING OF GENOMIC DNA SEGMENTS 


Sequences of each gene were obtained directly 
from PCR-amplified products. All oligonucleotide 
primers used for PCR and sequencing are listed in 
Table 2 along with their corresponding positions 
and orientations in the genes. The 5' portions of 
tub2, inclusive of introns IVSI, IVS2, and IVS3, 
were amplified using non-selective primers tub2- 
exonld-l and tub2-exon4u-1. A ca. 850 bp frag- 
ment of tef1, including the first five introns, was 
amplified using primers tefl-exonld-l and tefl- 
exon6u-1. A ca. 700 bp fragment of the act] gene, 
including introns 1—5, was amplified by PCR using 
primers actl-exonld-l and асі1-ехопби-1. Inter- 
nal regions of the PCR-amplified product were se- 
quenced with primers actl-exon3u-l and actl- 
exon3d-1. All amplification products were verified 
by agarose gel (0.8%) electrophoresis. Water blanks 
were included as controls to identify possible con- 
tamination. Successfully amplified products were 
purified for sequencing using Quiaquick spin col- 
umns (Quiagen, Chatsworth, California). PCR was 
performed as described above, except that the tem- 
perature program was as follows: one incubation at 
94°C for 9 min.; then 40 cycles of 94°C for 60 sec., 
55°С for 60 sec., and 72°C for 60 sec. Reaction 
mixtures were as described above. 

PCR products were sequenced by the Sanger 
method (Thomas & Kocher, 1993) with a Perkin- 
Elmer GeneAmp PCR System 2400 using rhoda- 


mine-labeled dideoxynucleotide triphosphates or 
BigDye Terminator Cycle sequencing kit (Perkin- 
Elmer). Sequence reaction products were analyzed 
using an ABI model 310 capillary electrophoresis 
genetic analyzer (Perkin-Elmer). Both DNA strands 
were sequenced. 


PHYLOGENETIC ANALYSIS 


Sequences were aligned using the PILEUP pro- 
gram in the Wisconsin Package (GCG, 1996). 
PILEUP parameters were adjusted empirically; a 
gap penalty of one and a gap extension penalty of 
zero resulted in accurate alignment of intron-exon 
junctions and in a reasonable alignment of the in- 
trons within the ingroup (Epichloé/Neotyphodium 
spp.) for individual analysis of tub2, act], and tefl 
genes, and for a combined data set using appended 
sequences of all three genes for representative iso- 
lates from each Epichloé species, as well as from 
outgroup taxa, Claviceps purpurea ATCC20102 and 
Echinodothis tuberiformis ATCC201937. Both out- 
group taxa are in the family Clavicipitaceae and are 
among the closest known relatives of Epichloé 
(Glenn et al., 1996). No ambiguous alignments 
were evident for the ingroups, but slight realign- 
ment of outgroup sequences were needed and were 
performed by eye. Maximum parsimony (MP) em- 
ployed the branch-and-bound option in PAUP* ver- 
sion 4.0 beta 4 (Swofford, 1998) for exact solutions. 
For parsimony analysis, character changes were un- 
weighted and unordered; gaps were treated as miss- 
ing information. Robustness of the data was esti- 
mated by bootstrap replications (Swofford et al., 
1996). All branches receiving 7096 or higher boot- 
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Table 3. 
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Linkage analysis between act], tub2, and tef1-alpha genes. Parental alleles represented by letters a (for 


E. typhina ATCC200736) and b (for E. typhina ATCC200740). Parental isolates were chosen based on differential 
polymorphism in each allele detectable by selective primers in PCR. 


Fungal isolate асі tefl 


E. typhina parental isolates used in cross 


ATCC200736 a a 
ATCC200740 b b 


Progeny isolates obtained 


E386.003 
E386.009 
Е386.026 
E386.027 
E386.042 
E386.043 
E386.045 
E386.057 
E386.084 
E386.180 


aa:bb:ab:ba 


c tW 


PPP PE POOH р 


Soe or im ue 


strap values were considered well supported. For 
distance-based analysis a number of different 
methods were employed to derive pairwise distanc- 
es. These included Jukes-Cantor and Kimura two- 
parameter distances [transition/transversion (ts/tv) 
= 2], and various gamma values to adjust for be- 
tween-site variation. Distance trees were inferred 
by neighbor-joining (NJ) implemented in PAUP* 
version 4.0 beta 4 (Swofford, 1998) with random 
taxon additions. All such distance trees were sim- 
ilar to the MP trees, and are therefore not shown. 
Maximum likelihood (ML) trees were constructed 
using PUZZLE in PAUP* version 4.0 beta 4. A 
Hasegawa-Kishino-Yano model (Swofford et al., 
1996) was employed with transition/transversion — 
2.0. PUZZLE uses quartet puzzling to create trees 
for all possible quartets of taxa via maximum like- 
lihood. These quartet trees serve as starting points 
to reconstruct a set of optimal trees for all taxa. The 
majority-rule consensus of these trees defines the 
quartet puzzling tree and shows well-supported 
branches (Strimmer & von Haeseler, 1996). Parti- 
tion homogeneity analysis was performed in PAUP* 
version 4.0 beta 4 on each two-gene combination 
in the combined sequence alignment. 


RESULTS 


LINKAGE ANALYSIS 


To check if the three genes selected for phylo- 
genetic analysis—namely, tub2, actl, and tef1— 
were unlinked, Mendelian segregation was assessed 
for variant alleles in an Epichloé typhina cross. Ten 


tub2 actl/tef1 actl/tub2 tef1/tub2 
a aa aa aa 
b bb bb bb 
b aa ab ab 
a ab aa ba 
a bb ba ba 
b bb bb bb 
b ab ab bb 
b ab ab bb 
a aa aa aa 
b ab ab bb 
a ab aa ba 
a ab aa ba 

2:2:6:0 4:1:4:1 1:4:1:4 


progeny were analyzed for alleles derived from the 
parents, ATCC200736 (male) and ATCC 200740 
(female). The occurrence of recombinant progeny 
was assessed for each gene pair. Ratios of parental 
to recombinant progeny genotypes indicated no 
close linkage between any of these genes in Ё. ty- 


phina (Table 3). 


MATING TESTS 


Past mating tests indicated that the currently de- 
scribed Epichloé species mostly correspond to in- 
terfertile groups (Leuchtmann & Schardl, 1998; 
Schardl & Leuchtmann, 1999). Exceptions were 
some instances of nonreciprocal interfertility among 
E. typhina isolates, and high fertility between Ё. 
clarkii and E. typhina. Several additional mating 
tests were conducted to establish species affinities 
of new isolates and mating tests that had not pre- 
viously been conducted for technical reasons. A 
new Epichloé isolate from Arrhenatherum elatius 
was fully compatible in matings with E. typhina 
stromata on Dactylis glomerata, as were the recip- 
rocal matings. The ascospore progeny were viable 
and exhibited segregation of parental allozymes. An 
isolate from Poa nemoralis was compatible with at 
least some strains from Brachypodium sylvaticum, 
but the matings exhibited low fertility with few asci 
containing ascospores and a very low germination 
rate of the ascospores. These matings were also 
confirmed by allozyme analysis. Epichloé typhina 
stromata from P. nemoralis mated successfully with 
Brachypodium pinnatum isolates, although fertility 
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was intermediate (moderate ascospore germination 
rate). Epichloé clarkii isolates from Holcus lanatus 
successfully mated with E. typhina stromata on Poa 
trivialis. Epichloé festucae fertilized with a new iso- 
late from Koeleria cristata reacted by forming pro- 
toperithecia, but no mature perithecia were ob- 
served and no ascospores were obtained. Epichloé 
festucae stromata on Festuca rubra mated success- 
fully with endophyte isolates from Lolium gigan- 
teum (= Festuca gigantea), Festuca valesiaca, Е 
glauca, and Е heterophylla. Matings were attempt- 
ed between stromata of E. amarillans isolates on 
Sphenopholis obtusata and spermatia of the E. 
amarillans ex-type (ATCC200744) from Agrostis 
hiemalis. Abundant perithecia gave rise to viable 
ascospores, but in the 10 ascospore isolates ana- 
lyzed, only maternal allozyme alleles were ob- 
served. 


PHYLOGENETIC RELATIONSHIPS BASED ON ACTI/ 
TUB2/TEF1 COMBINED DATA 


An MP tree (Fig. 2) was constructed from the 
combined data set of act], tub2, and tefl intron 
sequences, and included representatives of all de- 
scribed Epichloé species and two clavicipitaceous 
outgroups, Claviceps purpurea and Ech. tuberifor- 
mis. The tree shown in Figure 2 reflects the overall 
diversity within clades and described Epichloé spe- 
cies. However, results described later indicated 
some incongruences between gene trees due to re- 
ticulation involving E. typhina. For this reason the 
positions of E. typhina isolates did not reflect their 
positions in all gene trees. To test that the com- 
bined analysis gave a valid result for relationships 
between species, it was rerun with only two rep- 
resentatives of E. typhina, and again without any 
E. typhina sequences (data not shown). In all cases 
the interrelationships of the remaining taxa were 
the same. 

In bootstrap analysis (again using MP) of the 
combined data set most branches were strongly 
supported, but the placement of the root received 
only weak support (Fig. 3). This was not surprising 
because of the considerable sequence divergence 
between the ingroup and outgroup taxa, making un- 
ambiguous alignment difficult. Nevertheless, the in- 
ferred root was very close to the midpoint root of 
each individual gene tree (discussed later and 
shown in Figs. 4, 7, 11). Assuming its correct 
placement in the combined data analysis, the root 
split the Epichloé genus into two major groups (Fig. 
3): One group, comprising Epichloé amarillans, E. 
baconii, E. festucae, E. bromicola, E. elymi, E. gly- 
ceriae, and E. brachyelytri, is henceforth designated 
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E.amarillans 200743 So 
E.amarillans 200744 Agh 
E.baconii 76552 Ags 
E.baconii 90167 Agt 
E.festucae 90661 Рт 


















ww 20) chai 
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E.glyceriae 200747 Gs 
E.glyceriae 200755 Gs 
E.brachyelytri 200752 Bee 
E.brachyelytri 201561 Bee 
E.darkii 200742 НІ 
E.typhina 200736 Lp 
E.typhina 200851 Php 
E.typhina 201666 Ps 
E.typhina 200738 Ao 
E.typhina 200740 Dg 
E.sylvalica 200751 Bps 
E.typhina 200739 Врр 
E.sylvatica 200748 Bps 
E.typhina 201668 Pn 
E.typhina 201669 Pp 


жн 20102 
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201937 


Figure 2. МР tree generated from the combined act// 
tub2/tefl data set including representatives of each major 
Epichloé clade, Claviceps purpurea and Echinodothis tub- 
eriformis (outgroups). Branch-and-bound search resulted 
in 2 MP trees, differing slightly in the Е. festucae/E. ba- 
coni/E. amarillans clade. The alternative resolution of 
this clade is shown in the smaller box at right. The num- 
ber of parsimony-informative characters — 379; tree 
length — 1350 steps; consistency index — 0.8622; reten- 
tion index = 0.8502; rescaled consistency index = 
0.7331; outgroup rooting with basal polytomy. Bar repre- 
sents 20 inferred nucleotide substitutions. Letters follow- 
ing isolate number refer to abbreviated host designation. 
Full host names are listed in Table 1. 


the “main group” of Epichloé species. The second 
group, henceforth called the “E. typhina complex,” 
included E. typhina isolates from several host spe- 
cies, as well as E. clarkii (from Holcus lanatus) and 
E. sylvatica (from Brachypodium sylvaticum). A 
very strongly supported clade within the E. typhina 
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E.amarillans 200743 So 
E amarilans 200744 Agh 
E baconi 76552 Ags 
E baconi 90167 Ag 
E.festucae 90661 Рт 
Efestucae 028 F1 
E baconi 200745 Cv 
E bromicola 200749 Be 
100E bromicola 201558 Br 
X E bromicola 201559 Bb 
100 E chai 200850 Ev 
Eelyni 201551 Ec 
100 E gyceriae 200747 Gs 
E gyceriae 200755 Gs 
100 E brachyelytri 200752 Bee 
E brachyelyr 201561 Bee 
И Е clarki 200742 HI 
n E yphina 200736 Lp 
100E typhina 200851 Php 
E yphina 201666 Ps 
81 E.typhina 200738 Ao 
& E ypina 200740 Dg 
& pm syvatica 200751 Bos 
E ypina 200739 Bpp 
E sylvatica 200748 Bps 
и E.typhina 201668 Рп 
E yphina 201669 Pp 


Bootstrap on actt/tub2/teft 100 
combined data 


57 100 


main group of Epichloë spp. 


= 
crovvn clade 


E.typhina complex 


C.purpurea 20102 
Ech.tuberiformis 201937 


Figure 3. Bootstrap tree generated from the combined 
actl/tub2Aefl data set. Numbers at branches are the per- 
centage of trees containing the corresponding clade based 
on 1000 bootstrap replications. Only bootstrap values 
greater than 70% are shown and considered supportive of 
the clades indicated by the branches. 


complex was designated the crown clade (Fig. 3), 
and includes among others, E. typhina from Dac- 
tylis glomerata (which was among the original hosts 
on which Persoon described E. typhina and which 
was later selected as host of the lectotype) and E. 
clarkii from Holcus lanatus. 


B-TUBULIN GENE PHYLOGENY 


The tub2 data set included sequences from 59 
isolates of Epichloé spp. Phylogenies derived by MP 
(Figs. 4, 5), NJ (not shown), and ML quartet puz- 
zling (Fig. 6) were consistent with each other and 
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Eamarilans 200743 So 
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tub? MP 
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E.elymi 201554 ЕМ 
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E.elymi 201555 ЕМ 
E.elymi 201556 Eh * 
Eelymi 201551 Ec 
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Figure 4. Epichloë spp. tub2 gene tree based on max- 
imum parsimony (MP) analysis of introns 1—3. Shown is 
the single MP tree obtained by branch-and-bound search. 
Number of parsimony-informative characters — 82; un- 
informative characters — 21; tree length — 143 steps; 
consistency index = 0.8182; retention index = 0.9426; 
rescaled consistency index — 0.7712; midpoint root is at 
the left edge. Bar represents one inferred nucleotide sub- 
stitution. Asterisks (*) indicate the number of additional 
isolates (ranging from one to four), sampled from the same 
host genus, that had identical sequence. 


with the tub2 relationships previously reported 
among Epichloë species (Ѕсһагаї & Leuchtmann, 
1999). Both MP bootstrap analysis (Fig. 5) and ML 
quartet puzzling (Fig. 6) strongly supported the two 
major groups elucidated by the combined data set. 
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Figure 5. Bootstrap tree generated from Epichloé tub2 
gene sequence using MP (heuristic). Numbers at branches 
are the percentage of trees containing the corresponding 
clade based on 1000 bootstrap replications. Bootstrap val- 
ues greater than 70% are considered supportive of the 
clades indicated by the branches. For bootstrap analysis 
each set of identical sequences was represented twice. 
Asterisks (*) indicate the number of additional isolates 
(ranging from one to four), sampled from the same host 
genus, that had identical sequence. 


Members of most biological species tended to group 
together. For example, the isolates of E. amarillans, 
E. brachyelytri, E. festucae, E. bromicola, E. clarkii, 
E. elymi, and E. glyceriae grouped into exclusive 
clades. However, E. baconii from Agrostis spp. were 
not definitively grouped with Е. baconii from Cal- 
amagrostis villosa, the two E. sylvatica isolates did 
not definitively group in a monophyletic clade, and 
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Figure 6. Epichloë spp. tub2 gene tree based on max- 
imum likelihood (ML) PUZZLE analysis of introns 1—3. 
Tree is based on a Hasegawa-Kishino-Yano model with ts/ 
tv = 2 (kappa = 4.01068); number of puzzling steps = 
1000. Tree shown represents majority-rule consensus and 


tub2 puzzle 


numbers at branches are estimations of support for each 
branch, 


E. typhina exhibited apparent paraphyly with E. 
sylvatica and E. clarkii. 

The tub2 tree included sequences from several 
new isolates. An isolate from host species Koeleria 
cristata was placed in the Epichloé festucae clade, 
with perfect sequence identity to four other E. fes- 
tucae isolates. New isolates from Arrhenatherum 
elatius, Poa silvicola, and P. trivialis grouped to- 
gether with isolates from several other hosts in the 
crown clade of the E. typhina complex. New iso- 
lates from Poa pratensis and Brachypodium pin- 
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natum grouped with other isolates from these same 
hosts, in basal clades of the E. typhina complex. 


TRANSLATION ELONGATION FACTOR 1-0 GENE 
PHYLOGENY 


To assess the extent to which the tub2 gene tree 
reflected the evolutionary relationships among Ep- 
ichloé species and isolates, noncoding portions of 
two additional genes, tef] and act], were sequenced 
from most of the isolates. A portion of tef] contain- 
ing introns 1—4 was amplified and sequenced from 
representatives of the Epichloé species, and phy- 
logenetic trees were inferred by MP (Figs. 7, 8), NJ 
(not shown), and ML quartet puzzling (Fig. 10). The 
tefl MP trees resolved most of the same clades that 
were resolved in the tub2 trees, except that the po- 
sitions of E. glyceriae and E. brachyelytri were un- 
resolved in the tefl trees. Within the E. typhina 
complex tub2 and tef] phylogenies differed in the 
crown clade, as would be expected if clade mem- 
bers are interbreeding. In bootstrap analysis, all 
major clades that were strongly supported by the 
tub2 data were also strongly supported by tef1 (Fig. 
8). The tefl trees grouped isolates of each of the 
species, E. brachyelytri, E. bromicola, E. clarkii, E. 
elymi, E. festucae, and E. glyceriae. Furthermore, 
E. typhina isolates from Lolium perenne, Poa ne- 
moralis, and Brachypodium pinnatum grouped into 
exclusive clades, each with strong bootstrap sup- 
port. 


ACTIN GENE PHYLOGENY 


The act] MP trees were not as highly resolved 
as the other gene trees. This lack of resolution re- 
sulted in failure to successfully complete bootstrap 
analysis on the act] data set. Even so, all major 
clades indicated by tub2 and tefl analyses were 
supported in the act] consensus trees (Fig. 11). 
One conspicuous incongruence between the асі 
phylogeny and those constructed from both tub2 
and tefl was the placement of the isolates within 
the basal clades of the Epichloé typhina complex. 
Both tub2 and tef] analysis placed the E. typhina 
isolates from Poa nemoralis and P. pratensis togeth- 
er in the most basal clade, whereas the act] trees 
placed the E. typhina isolates from Brachypodium 
pinnatum together with one of the Е. sylvatica iso- 
lates (ATCC 200751) from Bp. sylvaticum as the 
most basal clade (Figs. 11, 12). Such incongruenc- 
es between gene trees in placement of taxa sug- 
gested possible genetic reticulation due to inter- 
breeding within the E. typhina complex. 
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Figure 7. Epichloë spp. tefl gene tree based on МР 
analysis of introns 1—4. Branch-and-bound search resulted 
in 72 MP trees. Number of parsimony-informative char- 
acters = 110; tree length = 196 steps; consistency index 
= 0.8571; retention index = 0.9506; rescaled consistency 
index = 0.8148; midpoint root is at the left edge. Bar 
represents five inferred nucleotide substitutions. Asterisks 
(*) indicate the number of additional isolates (ranging 
from one to two), sampled from the same host genus, that 
had identical sequence. Some sequences separated on 
zero-length branches differ only in indels, which were 
treated as missing information in this analysis. 


TESTS FOR TREE INCONGRUENCES 


To circumscribe interbreeding populations we in- 
vestigated reticulation within clades. Interbreeding 
should lead to incongruences between gene trees 
in the placement of isolates within the interbreed- 
ing population. Most clades corresponding to inter- 
fertile groups had too little sequence variation to 
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Figure 8. Epichloé spp. tefl gene tree based on MP 
analysis of introns 1—4. Branch-and-bound search resulted 
in 72 MP trees. Tree shown was chosen to represent var- 
iation in branching orders among the MP trees. The most 
important differences are the positions of E. brachyelytri 
and E. glyceriae (boxed). Number of parsimony-informa- 
tive characters = 110; tree length = 
tency index — 0.8571; retention index — 0.9506: rescaled 
consistency index — 0.8148; midpoint root is at the left 
edge. Bar represents five inferred nucleotide substitutions. 
Asterisks (*) indicate the number of additional isolates 
(ranging from one to two), sampled from the same host 
genus, that had identical sequence. Some sequences sep- 
arated on zero-length branches differ only in indels, which 
were treated as missing information in this analysis. 


196 steps; consis- 


make such an assessment. However, three such 
groups were highly diverse in sequences: Epichloé 
baconii, E. sylvatica, and the Е. typhina/E. clarkii 
group. In all three gene trees (Figs. 4, 7, 11) the 
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Figure 9. Strict consensus cladogram from the 72 MP 
trees derived from a branch-and-bound search of tef] se- 
quence. Numbers at branches are the percentage of trees 
containing the corresponding clade based on 500 boot- 
strap replications. Bootstrap values greater than 70% are 


considered supportive of the clades indicated by the 
branches. For bootstrap analysis each set of identical se- 
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quences was represented only once. 


E. baconii isolate from Calamagrostis villosa was 
well separated from those with Agrostis spp. hosts. 
Thus, there was no conflict between the gene trees 
that would suggest intermating between Е. baconii 
in different clades. There were, however, incongru- 
ences between gene trees in the E. typhina com- 
plex. Specifically, there were differences in the 
crown clade between all gene trees, and there were 
differences in basal clades of the complex when the 
act] tree was compared to the tub2 and tefl trees. 
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Figure 10. Epichloë spp. tefl gene tree based on ML 
puzzle analysis of introns 1—5. Tree is based on a Hase- 
gawa-Kishino-Yano model with ts/tv = 2 (kappa = 
3.97409); number of puzzling steps = 1000. Tree shown 
represents majority-rule consensus, and numbers at 
branches are estimations of support for each branch. 


lf puzzle 


These incongruences and the involvement of these 
clades were investigated by partition homogeneity 
tests. 

The partition homogeneity test compares phylo- 
genetic information from two genes to determine 
whether there is statistically significant evidence of 
incongruence between sequence data sets (Р = 
0.05). Since we were interested in reticulation (pos- 
sible intermating) within the genus Epichloé, we re- 
moved the outgroup taxa (C. purpurea and Ech. tub- 
eriformis) from the analysis. Comparison of act] 
with tef] resulted in significant incongruence (P = 
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Figure 11. Strict consensus cladogram from the 1946 
MP trees derived from a branch-and-bound search of act] 
sequence. Number of parsimony-informative characters = 
80; tree length = 149 steps; consistency index = 0.8792; 
retention index = 0.9430; rescaled consistency index = 
0.8291; midpoint root is at the left edge. Asterisks (*) 
indicate the number of additional isolates (ranging from 
one to three), sampled from the same host genus, that had 
identical sequence. 
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0.03). To determine which taxa may be reticulated, 
thus the source of incongruence, we removed cer- 
tain members of the E. typhina complex and re- 
analyzed the data. When all of the crown clade was 
removed except E. typhina from Dactylis glomerata 
and E. clarkii from Holcus lanatus, and the test was 
rerun, there was no longer significant incongruence 
between act] and tef] trees (P = 0.20). Similarly, 
removal of all E. typhina isolates except one from 
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Figure 12. Epichloë spp. act] gene tree based on ML 
puzzle analysis of introns 1—4. Tree is based on a Hase- 
gawa-Kishino-Yano model with ts/ty = 2 (kappa = 
3.98002); number of puzzling steps = 1000. Tree shown 
represents majority-rule consensus and numbers at 
branches are estimations of support for each branch. 


D. glomerata, while leaving Е. clarkii and Е. syl- 
vatica, gave a nonsignificant test result (P — 0.47). 

Partition homogeneity analysis of tub2 and tefl 
gave similar results. When the test included all Ep- 
ichloé taxa, there was a significant incongruence 
between these two gene sets (P — 0.01). Removal 
of all crown clade isolates except E. typhina from 
Dactylis glomerata and E. clarkii from Holcus lan- 
atus resulted in nonsignificance (P — 0.71). Elim- 
ination of all E. typhina isolates except the D. 
glomerata isolate also resulted in nonsignificance 


(P — 0.59). 
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Inclusion of all Epichloé isolates in the act1/tub2 
test indicated incongruence (P — 0.01). In contrast 
to the other gene comparisons, when the test was 
run without the crown clade E. typhina isolates it 
still indicated incongruence (P — 0.02). However, 
if taxa from Poa pratensis, P. nemoralis, and Bra- 
chypodium pinnatum were also pruned from the 
test, there was no significant incongruence (P — 
0.23). Thus, evidence of reticulation in this last 
comparison was only lost once Е. typhina isolates 
from both the crown and basal clades were re- 
moved, suggesting reticulation throughout the E. ty- 
phina complex. 


DISCUSSION 


PHYLOGENETIC AND BIOLOGICAL SPECIES 


The main focus of this study was to examine the 
relationship between biological and phylogenetic 
species, as well as the role of ecological speciali- 
zation in genetic isolation of populations of Epich- 
loé species. To address these issues, we sought ad- 
ditional genes to compare to the tub2 gene trees 
previously inferred (Schardl et al., 1997). Since ge- 
netic linkage tests indicated that tefl, act1, and 
tub2 were not closely linked, at least in E. typhina, 
these were considered appropriate for gene tree 
comparisons. Taxonomic sampling was comparable 
in scope to our recent tub2 phylogenetic analysis 
(Schardl & Leuchtmann, 1999), with all 10 de- 
scribed species of Epichloé sensu stricto represent- 
ed by multiple isolates. The trees inferred from tefl 
and act] intron sequences largely corroborated the 
relationships within the genus Epichloé previously 
elucidated by tub2 phylogeny (Tsai et al., 1994; 
Schardl & Leuchtmann, 1999). Analysis of com- 
bined data with outgroups also corroborated, with 
moderate support, the midpoint root inferred for 
each of the three gene trees. The tree thus rooted 
divides the genus into two major clades. One clade 
comprises most of the currently accepted Epichloé 
species, three of which (E. baconii, E. bromicola, 
and Е. festucae) are associated with Eurasian grass 
hosts (e.g., Festuca spp. and Bromus spp.), and four 
(E. amarillans, E. brachyelytri, E. elymi, and Е. 
glyceriae) with North American grasses (e.g., Ely- 
mus spp. and Agrostis spp.). This major Epichloé 
clade we refer to as the main group. The other ma- 
jor clade, here termed the Е. typhina complex, was 
associated with a broad range of Eurasian grasses 
and included the described species E. typhina, E. 
clarkii, and E. sylvatica. For the most part, the main 
group of Epichloé comprised clear-cut phylogenetic 
species corresponding to established interfertility 
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groups and morphologically based species. A pos- 
sible exception, E. baconii, is discussed later, as is 
the possibility of cryptic biological species within 
morphospecies E. festucae and E. amarillans. How- 
ever, the E. typhina complex poses the most diffi- 
culty in reconciling phylogenetic and biological 
species concepts. 

Biological species can be defined as potentially 
interbreeding populations that are reproductively 
isolated from other mating populations (Claridge et 
al., 1997; Mayr, 1940). Recent descriptions of Ep- 
ichloé species have considered the results of mating 
tests to identify potentially interbreeding groups 
(Leuchtmann et al., 1994; Leuchtmann & Schardl, 
1998; Schardl & Leuchtmann, 1999; White, 1993, 
1994). However, such tests do not necessarily in- 
dicate actual interbreeding populations. This is a 
particular concern if parents are so specialized to 
different ecological niches, such as hosts, as to pre- 
clude survival of their progeny. A phylogenetic ap- 
proach to identify effectively interbreeding popu- 
lations has recently been applied to fungi for which 
the sexual cycle is unknown (Geiser et al., 1998; 
Burt et al., 1996). The approach is to identify the 
smallest well-defined phylogenetic units in com- 
bined analysis of multiple genes. Such clades are 
suggested to indicate reproductively isolated pop- 
ulations (Geiser et al., 1998). One assumption is 
that phylogenetic species, so defined, are equiva- 
lent to biological species. In our study we have 
looked only at known sexual Epichloé species and, 
furthermore, have identified the interfertility groups 
in the genus. This allows a check on the relation- 
ships between phylogenetic and biological species. 

We looked for reticulation that would give rise to 
disparities between tree topologies inferred from 
different, unlinked genes. Reticulation indicates 
that either the population has a history of recom- 
bination characteristic of an interbreeding gene 
pool or that much of the extant polymorphism was 
inherent in the clade’s ancestral gene pool. We 
point out that this analysis cannot indicate how re- 
combination occurs, how often it occurs, or when it 
last occurred. Genetic tests indicated that the tefl, 
actl, and tub2 genes independently segregated in 
Epichloé typhina, so we considered these genes ap- 
propriate for the analysis. Tests of linkage were 
most feasible in E. typhina because this species 
had sufficient sequence diversity of all three genes 
for the appropriate primers to be developed and 
PCR tests to be conducted. An assumption was that 
these genes are likewise not closely linked in other 
Epichloé clades. We consider this a reasonable as- 
sumption because, since they are not linked in E. 
typhina, it is unlikely that all three would be so 
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closely linked in other clades as to affect inferences 
about genetic reticulation. More importantly, the 
most apparent conflicts between phylogenetic and 
biological species were associated with the E. ty- 
phina complex. 

Likely reticulation was identified by conducting 
partition homogeneity tests both with all taxa in- 
cluded and with most members of potentially retic- 
ulating clades removed. Since all three gene trees 
were obviously concordant for the main Epichloé 
group, we focused on the E. typhina complex, 
where the different gene trees gave different topol- 
ogies. Not surprisingly, when all taxa were included 
the test indicated significant conflict between the 
trees. In two comparisons (act1/tef1 and tub2/tef1) 
significance was lost when the Е. typhina crown 
clade (including E. clarkii) was reduced to two rep- 
resentatives. This result strongly suggests that the 
crown clade comprises or is a component of a re- 
ticulated gene pool. Only when the F. typhina com- 
plex was further pruned by removal of most basal 
taxa did the test indicate congruent act] and tub2 
trees. This suggests that there was also reticulation 
in the basal taxa, and inspection of the gene trees 
indicates that the positions of E. sylvatica isolates 
differ in the different trees. The positions of Е. syl- 
vatica isolates in the three gene trees were a par- 
ticular surprise, as was the lack of fertility between 
E. sylvatica and Brachypodium pinnatum-associat- 
ed E. typhina despite their close relationships 
(Leuchtmann & Schardl, 1998). It would appear 
that Е. sylvatica has emerged as a distinct species 
from within Ё. typhina, or that E. sylvatica has re- 
combined with Е. typhina in the course of its evo- 
lution prior to, or despite, emergence of intersteril- 
ity barriers. 

It is possible that this conflict between biological 
and phylogenetic species may be a problem of lin- 
eage sorting in the Epichloé typhina complex. Lin- 
eage sorting arises due to polymorphisms present 
in an ancestral mating population from which dif- 
ferent alleles are fixed in extant lineages, and can 
result in incorrect phylogenetic inferences being 
drawn as to the true relationships between extant 
lineages (Lyons-Weiler & Milinkovitch, 1997). An- 
other possibility is that there has been parasexual 
recombination involving Ё. typhina and Е. sylvatica 
(Caten, 1981). This is plausible considering the or- 
igins of many asexual endophytes as hybrids of 
multiple Epichloé species. However, the hybrid 
asexual species generally have large portions of the 
genomes of their ancestors, thus containing multi- 
ple copies of most genes (Kuldau et al., 1999). 
There is no indication of this in Epichloé species 
from analysis of tub2, tef1, and act1, though occa- 
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sional multiband phenotypes have been observed 
in isozyme studies of E. typhina (Leuchtmann & 
Schardl, 1998). Even if parasexual recombination 
occurs in the Е. typhina complex, the close rela- 
tionship between E. sylvatica ATCC 200751 and E. 
typhina from Brachypodium pinnatum is unexpect- 
ed in that all three gene trees give this same rela- 
tionship, and isozyme analysis gives a similar result 
(Leuchtmann & Schardl, 1998). In stark contrast, 
E. typhina from Bp. pinnatum did not share any 
identical or highly similar gene sequences with Е. 
typhina isolates from any other host, as would have 
been expected from either sexual or parasexual re- 
combination. Although we cannot exclude either 
lineage sorting or parasexual recombination to ex- 
plain the apparent paraphyly of E. typhina to Е. 
sylvatica, as well as polyphyly of both biological 
species, these relationships nevertheless amount to 
a conflict between the phylogenetic and biological 
species concepts. 

With the exception of the E. typhina complex, 
most phylogenetic species of Epichloé had little ge- 


netic diversity, and there is little or no evidence of 


gene flow (reticulation) between phylogenetic spe- 
cies. Previous mating test results (Schardl & 
Leuchtmann, 1999) have shown that many phylo- 
genetic species in the main Epichloé group also 
conform to a biological species concept; there was 
intermating only between isolates from the same or 
closely related hosts, often within one host genus. 
Yet, even in the main group there were some cases 
where phylogenetic species did not precisely pre- 
dict interfertility groups or vice versa. Epichloé ba- 
conii isolates from host genera Agrostis and Cala- 
magrostis, though capable of intermating, never 
definitively grouped in a monophyletic clade in any 
gene phylogram (this was also the observation 
based on allozyme cluster analysis; Leuchtmann & 
Schardl, 1998). However, it is possible that when- 
ever Agrostis- and Calamagrostis-associated popu- 
lations intermate in nature their progeny are lost 
due to lack of a compatible host. Such a situation 
would isolate the populations into different biolog- 
ical species despite their observed interfertility in 
experimental matings. More sampling of the Е. ba- 
conii populations is needed to assess if they actu- 
ally represent distinct biological species or, alter- 
natively, if this is another conflict between 
phylogenetic and biological species concepts. 

The converse situation apparently arises with 
Epichloé amarillans and E. festucae. An E. amar- 
illans isolate on Sphenopholis obtusata had limited 
or no interfertility with an isolate from Agrostis hie- 
malis. Although mating proceeded as expected, 
progeny showed no evidence of genetic contribution 
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from the male. Though this does not preclude lim- 
ited genetic exchange, the result suggests that these 
host-associated populations have a degree of ge- 
netic isolation. Since phylogenetic analysis placed 
these populations into distinct clades they may in- 
deed represent different biological species. Like- 
wise, E. festucae from Festuca rubra failed to mate 
with an isolate from Koeleria cristata even though 
they had identical gene sequences. This test must 
be considered preliminary pending reciprocal mat- 
ing tests. Nevertheless, the result points up a po- 
tential limitation that if sequence evolution does 
not adequately keep pace with biological speciation 
it may be difficult to identify species on a strictly 
phylogenetic basis. 


GENE FLOW WITHIN A BIOLOGICAL SPECIES 


The relationships among host-associated popu- 
lations in the Epichloé typhina complex is intrigu- 
ing. In keeping with previous results from isozyme 
analysis (Leuchtmann & Schardl, 1998) we found 
that certain host grasses harbor genetically diverse 
endophyte populations, whereas other hosts harbor 
populations of genetically similar individuals. 
Among the more diverse are the E. typhina popu- 
lations in Dactylis glomerata and Phleum pratense, 
as well as E. sylvatica from Brachypodium sylvati- 
cum. In contrast, isolates from other host-associated 
fungal populations had similar or identical allozyme 
profiles and DNA sequences. Little or no sequence 
variation was observed within Е. typhina popula- 
tions in many of its grass hosts: Anthoxanthum 
odoratum (based on tub2), Brachypodium pinnatum 
(tub2 and tefl), Lolium perenne (tub2, tefl, and 
асі1), Poa nemoralis (tub2, tefl, and act1), and P. 
pratensis (tub2); as well as among E. clarkii from 
Holcus lanatus (tub2 and tef1). Since most of these 
host-associated populations are currently repre- 
sented by a small sample size it is unclear whether 
they might be substantially more diverse than so 
far indicated. Still, the genetic uniformity of Е. 
clarkii seems well substantiated because this spe- 
cies Is represented by five isolates from Britain and 
Switzerland, and both mating types are included. 
Even though most F. clarkii isolates are readily in- 
terfertile with E. typhina on several other hosts, E. 
clarkii seems to maintain genetic and even mor- 
phological distinction from the rest of the E. typhi- 
na complex (White, 1993; Leuchtmann & Schardl, 
1998). However, the position of E. clarkii relative 
to the various Ё. typhina isolates is not consistent 
among gene trees, so that E. clarkii does not fit the 
criterion for a phylogenetic species. It is possible 
that E. typhina on other grasses—such as Dactylis 
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glomerata—actually mates with E. clarkii in na- 
ture, but that the progeny do not colonize or survive 
on Holcus lanatus, host grass for E. clarkii. Perhaps 
only E. clarkii X E. clarkii progeny can colonize 
this host, leading to extreme inbreeding but not 
necessarily genetic isolation. Thus, E. clarkii may 
contribute to the gene pool of certain E. typhina 
populations, but Е. typhina rarely or never contrib- 
utes to the E. clarkii (i.e., Н. lanatus-associated) 
gene pool. The other low-diversity host-associated 
populations may have similar relationships to the 
more diverse ones. 

Is incipient speciation occurring in the Epichloé 
typhina complex? One might expect speciation 
based on host-associated populations given the 
strong association of other Epichloé species with 
individual host genera or tribes. Epichloe sylvatica 
appears to be an example, since it is restricted to 
one host, Brachypodium sylvaticum. Yet even if 
such species emerge, it appears that they will not 
correspond to such clear-cut phylogenetic clades as 
do the other Epichloé species. This, also, is illus- 
trated by the apparent polyphyly of E. sylvatica. 
Overall, no phylogenetic species obviously emerge 
from host-associated populations in the Ё. typhina 
complex. However, it appears that the crown clade 
of E. typhina/E. clarkii may have effectively spe- 
ciated from the more basal clades. All indications 
of reticulation were either within the crown clade 
or among the basal clades, but not between the 
crown and basal clades. This suggests that the 
crown clade is an incipient species that has not 
fully established genetic barriers to mating outside 
the clade (many matings between basal and crown 
clade members are fertile; Leuchtmann & Schardl, 
1998). Furthermore, even though one genus (Poa) 
is represented among hosts of both basal and crown 
clades, the different clades consistently comprise 
isolates from different host species. Epichloé spe- 
cies in the crown clade infect P. silvicola and Р. 
trivialis, whereas members of the more basal clades 
infect P. pratensis and P. nemoralis. 


PREVAILING TRANSMISSION MODES AND EVOLUTION 


The apparently different evolutionary processes 
between the Epichloé typhina complex and the 
“main group" containing the other seven Epichloé 
species may relate to their different host interac- 
tions, namely whether they are more benign and 
vertically transmissible or more virulent (antago- 
nistic) and reliant on horizontal transmission. Most 
members of the E. typhina complex are antagonistic 
(strong chokers) in that they shut down almost all 
host seed production. They must therefore rely on 
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horizontal transmission. As discussed above, spe- 
ciation in this complex is not clear cut. In contrast, 
most of the species in the main group balance ver- 
tical and horizontal transmission, and, as partial 
chokers, have more benign or even mutualistic ef- 
fects on host plants. In this group speciation follows 
a cladistic stereotype and most species boundaries 
are unambiguous. 

Species in the main Epichloé group with sub- 
stantial seed transmission include E. brachyelytri, 
E. amarillans, E. festucae, and E. elymi (Leucht- 
mann & Schardl, 1998; Leuchtmann et al., 1994; 
Schardl & Leuchtmann, 1999). The situation for E. 
bromicola is unique in that it is a strong choker on 
Bromus erectus, but is strictly seed transmitted on 
other Bromus species (Leuchtmann & Schardl, 
1998). It is possible that E. bromicola has both 
transmission modes on an unidentified host. Only 
two main group species, Ё. glyceriae and Е. bacon- 
ii, are strong chokers for which vertical transmis- 
sion is unknown. The deep rooting of E. brachye- 
lytri (based on combined sequence analysis; Fig. 3) 
and the preponderance of vertically transmitted 
species suggest that this is an ancestral trait of the 
main Epichloé group. We previously proposed that 
the loss of seed transmissibility in E. glyceriae 
might have been associated with a host jump, since 
without this fungal species the relationships among 
the main Epichloë group largely reflect host rela- 
tionships (Schardl et al., 1997). However, since the 
relevant branches were unresolved in tef] and act] 
trees, it remains possible that E. glyceriae and the 
rest of the main group have tracked the evolution- 
ary paths of their hosts. 

Even with inconsistent phylogenetic tracking it 
is likely that host specialization played a role in 
Epichloé speciation. Since vertical transmission 
provides a dependable route for dissemination 
there is little benefit expected for broadening host 
range. This is particularly true if there is a tradeoff 
between broad host range and the balance of ver- 
tical and horizontal transmission in these fungi. We 
have suggested that such a balance of host and fun- 
gal development is unlikely to survive host jumps, 
and that this is why the more mutualistic Epichloé 
species seem to track host phylogeny (Schardl et 
al., 1997). Here we suggest that this is also why 
there are clear boundaries between host-associated 
species in the main Epichloé group. As further ev- 
idence that benign or mutualistic interactions have 
promoted host-based speciation, the phylogenetic 
pattern is substantially different for the more an- 
tagonistic species, especially the E. typhina com- 
plex. Models of host-parasite interactions predict 
that high virulence levels (antagonism) should cor- 
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relate with rapid contagious spread (Yamamura, 
1993). No comparative study has been done on 
transmission rates of Epichloé species. However, 
horizontal transmission of E. typhina appears to be 
very rapid in Dactylis glomerata populations, as in- 
dicated by a recent epidemic in Oregon (Pfender 
& Alderman, 1999). We speculate that high viru- 
lence and reliance on horizontal transmission se- 
lects for broad host range and limits speciation in 
the E. typhina complex. There are two reasons why 
E. typhina may be particularly subject to this kind 
of selection. First, since this species depends 
heavily on horizontal transmission the availability 
of hosts within a short range of a parent stroma 
(specifically, the dispersal area for ascospores) is 
important for its propagation. Second, it can partic- 
ipate in reducing the host population, and indeed 
it will be a factor in the reduction of available hosts 
by reducing flowering. Whenever local host popu- 
lations decline either by effects of E. typhina or 
other factors, there would be increased selection on 
the fungus to colonize other hosts. Such a selective 
effect is not expected for more benign Epichloé spe- 
cies because their assured vertical dissemination 
obviates any immediate reliance on horizontal 
transmission to neighboring plants. 


ORIGINS OF THE MUTUALISTIC ENDOPHYTES 


The main group of Epichloé species represents 
the emergence of a mechanism to maximize mutual 
benefit to host and symbiont. This mechanism is a 
balanced interaction whereby individual plant-fun- 
gus symbiota can simultaneously exhibit sexual 
propagation of both partners. Although the neces- 
sary components of such interactions are repre- 
sented among related Clavicipitaceae (Clay & 
Jones, 1984; Kover & Clay, 1998), it is only in the 
main Epichloé group that such a balance is consis- 
tently evident (Schardl & Clay, 1997). In contrast, 
the more antagonistic effects of the E. typhina com- 
plex are similar to those of many other plant-as- 
sociated Clavicipitaceae, particularly Balansia spe- 
(Diehl, 1950; Kuldau et al., 1997). 
Furthermore, most of the E. typhina complex, like 
most Balansia species, fail to transmit in host 
seeds. Although Ё. typhina on Роа nemoralis can 
transmit vertically, it tends to choke the vast ma- 
jority of tillers and thereby severely limit both host 
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seed production and its potential for vertical trans- 
mission. Likewise, although Е. sylvatica was pre- 
viously considered to have balanced host interac- 
tions, more detailed examination has indicated that 
stroma-forming genotypes are actually strong chok- 
ers. What were previously thought to be balanced 
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interactions were simultaneous infections, whereby 
individual Brachypodium sylvaticum plants had two 
endophytes, one sexual (E. sylvatica) and the other 
asexual (Meijer & Leuchtmann, 1999). Though 
closely related to E. sylvatica, these asexual endo- 
phytes represent the opposite extreme in that they 
have sacrificed all capacity for horizontal spread 
and are strictly vertically transmitted. Although this 
may well benefit the host, it is nevertheless appar- 
ent that the balanced interactions characterizing 
much of the Epichloé main group species are absent 
in the E. typhina complex as they are among other 
Clavicipitaceae. 

Symbioses of asexual endophytes with grasses 
represent a third class of fungal/host interactions 
distinct from the balanced or antagonistic associa- 
tions involving sexual Epichloé species. Many 
grasses in subfamily Pooideae possess asexual de- 
rivatives of Epichloé species, which, by current tax- 
onomic convention, are classified in a distinct ge- 
nus in the Fungi Imperfecti. Several such species 
have been formally described and placed in genus 
Neotyphodium (Glenn et al., 1996). These endo- 
phytes are often beneficial to their hosts and cause 
no disease, so they are commonly regarded as mu- 
tualists. We are currently analyzing tub2, tef1, and 
actl noncoding regions for several Neotyphodium 
species, some formally described and others re- 
maining to be described. Results from this tub2 
analysis indicate that several have hybrid origin 
(Schardl et al., 1994; Tsai et al., 1994; Moon et al., 
2000), a conclusion also supported by analysis of 
tefl and act] genes (K.D.C. & C.L.S., unpublished 
obs.). As we expand the number of genes and spe- 
cies analyzed it is becoming apparent that hybrid 
origins are very common among Neotyphodium 
spp., though not universal. For example, the asex- 
ual N. coenophialum, a common endophyte of Lol- 
ium arundinaceum (— Festuca arundinacea), has 
three Epichloé species in its pedigree (Tsai et al., 
1994); whereas N. lolii, the most common endo- 
phyte of Lolium perenne, seems to be a direct de- 
scendant of E. festucae. Those Epichloë species in- 
the 
phylogenetic diversity of the genus, including E. 
typhina, E. bromicola, E. elymi, E. amarillans, E. 
festucae, and a relative of E. baconii. The hybrid- 


volved in known hybridizations span 


izations have led to diploidy, polyploidy, or exten- 
sive aneuploidy in most Neotyphodium species 
studied to date, and for this and other reasons are 
likely somatic events rather than hybrid matings 
(Schardl et al., 1994; Kuldau et al., 1999). Thus, 
traditional species concepts are violated, albeit in 
different ways, both by endophytes that only trans- 
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mit vertically and those that only transmit horizon- 
tally. 


CONCLUSIONS 


The most striking result from our phylogenetic 
studies of Epichloé species and their asexual rela- 
tives is that different evolutionary patterns are as- 
sociated with different modes of transmission to 
new host plants. A cladistic pattern, usually as- 
sumed to describe speciation in most organisms, 
was only apparent for one of the two major lineages 
in the Epichloé phylogeny. Surprisingly, this was the 
group that had an exquisite balance of horizontal 
and vertical transmission uncharacteristic of the 
fungal family and, indeed, uncharacteristic of most 
pathogens and symbionts. In contrast, cladistic evo- 
lution is less evident both for strictly vertically 
transmitted and strictly horizontally transmitted 
species. Interspecific hybridization is often evident 
in the origins of the vertically transmitted, asexual 
symbionts (Schardl et al., 1994; Tsai et al., 1994). 
Further, the E. typhina complex—dominated by 
strains that only transmit horizontally—lacks clear- 
cut genetically isolated groups (this study). We 
speculate that high virulence and dependence on 
horizontal transmission might counteract host spe- 
cialization and, in turn, promote intraspecies di- 
versity and delay speciation. Delayed speciation 
should amplify lineage sorting problems that cause 
conflict between the biological and phylogenetic 
species concepts, as discussed above. If a reliance 
on horizontal transmission delays speciation as we 
speculate, an obvious prediction is that other highly 
virulent pathogens will exhibit conflicts between bi- 
ological and phylogenetic species. To our knowl- 
edge, this has not yet been tested. A good relation- 
ship between these species concepts was noted for 
the Gibberella fujikuroi group (O'Donnell et al., 
1998), but these fungi also exhibit both horizontal 
and vertical transmission (Bacon & Hinton, 1996). 
Detailed phylogenies and population surveys are 
clearly needed for more virulent pathogens, organ- 
isms that cause the most problems in human health 
and the food supply. 
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